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DECOMPOSING  THE  JOINT  DISTRIBUTION  OF  SEVERAL 


LINEAR  COMBINATIONS  OF  UNIFORM  SPACINGS 
By 

Fred  Huffer 


1.  The  Main  Result. 

Place  n-1  independent  uniformly  distributed  points  on  the  unit 
interval.  These  points  divide  the  unit  interval  into  n  pieces  whose 
lengths  (called  spacings)  are  denoted  S^.S^,...^  .  The  random  vector 
S  *  (S^»S2»*« • »S  )’  is  uniformly  distributed  on  the  simplex  determined 
by  S^+S^ +  • * •  +  *  1.  It  follows  that  the  are  exchangeable  random 

variables.  For  further  information  on  spacings  see  Pyke  (1965). 

Let  T  be  a  k*n  real  matrix.  We  wish  to  investigate  the  distri¬ 
bution  of  the  random  vector  TS.  Let  P (T)  denote  the  probability 
measure  of  TS  so  that  Prob{rSeA)  =  (P(T))(A).  We  shall  present  an 
identity  which  in  some  cases  allows  us  to  express  P(T)  as  a  sum  of 
simpler  terms. 

For  any  C  ,  define  r.  _  to  be  the  k*n  matrix  obtained  by 

*  »  4 

replacing  the  i—  column  of  T  by  £  .  Let  e  eRn  be  the  vector  of 
ones,  e'  ■  (1,1,..., 1).  Our  main  result  is  the  following. 

Theorem:  Suppose  d'  =  (d„ ,d„, . . . ,d  )  satisfies  e'd  ■  1.  Let 
1 "  ’ "  i  z  n 

£  *  Td.  Then 

n 

P(D  -  l  d  p(r  ). 
i-1 


1 


This  theorem  is  proved  in  section  2.  Section  3  contains  many  examples 
using  the  theorem  to  find  probabilities  in  situations  involving  one  or  sev¬ 
eral  linear  combinations  of  spacings.  Other  approaches  to  these  problems 
have  been  given  in  the  literature.  A  number  of  authors  have  obtained  the 
distribution  of  a  single  linear  combination  of  spacings  which  we  deal  with 
in  example  1.  Finding  this  distribution  is  equivalent  to  determining  the 
volume  of  the  intersection  of  a  simplex  with  a  half-space.  Dempster  and 
Kleyle  (1968)  employ  an  inclusion-exclusion  argument  to  decompose  this 
region  into  simplices  and  thus  obtain  the  volume  in  the  case  when  the  coef¬ 
ficients  of  the  linear  combination  are  all  different.  Weisberg  (1971) 
generalized  their  result  to  arbitrarily  coincident  coefficients.  Varsi 
(1973)  gave  an  algorithm  for  computing  this  volume  based  on  a  different 
method  of  dissecting  the  region  into  simplices.  Ali  (1973)  obtains  the 
distribution  by  inversion  of  its  characteristic  function.  A  more  elementary 
analytic  argument  (avoiding  inversion  of  transforms)  is  given  by  Gerber 
(1981).  Other  closely  related  results  concern  the  distribution  of  a  linear 
combination  of  order  statistics  from  the  uniform  distribution  (Ali  (1969)) 
and  the  distribution  of  a  linear  combination  of  independent  exponential 
random  variables  (Ali  and  Obaidullah  (1982),  Mathai  (1983)).  The  papers  by 
Ali  (1973),  Gerber  (1981)  and  Ali  and  Obaidullah  (1982)  all  make  use  of 
properties  of  divided  differences  as  does  the  present  paper. 

The  problem  of  finding  the  joint  distribution  of  several  linear  com¬ 
binations  of  uniform  spacings  has  received  less  attention.  Some  special 
cases  were  dealt  with  by  Kleyle  (1967);  see  the  remarks  in  section  3  of 
Dempster  and  Kleyle  (1968)  and  Theorem  A. 2  of  Kleyle  (1971).  The  closely 


related  problem  of  finding  the  joint  distribution  of  linear  combinations 
of  order  statistics  from  the  uniform  distribution  was  considered  by  Ali 
and  Mead  (1969)  who  obtained  results  by  inverting  the  joint  characteristic 
function. 

There  is  a  large  body  of  work  dealing  with  probabilities  of  clustering 
on  the  line  or  circle,  probabilities  of  multiple  coverage  of  the  line  or 
circle,  and  the  distribution  of  the  scan  statistic  and  other  tests  for 
uniformity.  Specialized  methods  have  been  developed  to  calculate  these 
probabilities.  See  Naus  (1982)  for  references.  These  problems  frequently 
involve  calculating  P{X^  >  t  for  1  s  i  s  k)  or  P{X^  <  t  for  lsisk)  where 
(Xj^, . . .  ,X^)  "  =  TS  .  Here  r  is  some  matrix  whose  entries  are  all  zeros 
and  ones  which  are  usually  arranged  according  to  a  simple  pattern.  For 
example,  Glaz  and  Naus  (1983)  find  the  variance  of  the  number  of  clusters 
on  the  line  by  calculating 

k  j+k-1 

(*)  P<  I  S.<t,  l  S  <t> 

i=l  1  i=j  1 

for  arbitrary  j,  k  and  t.  Examples  3  through  6  illustrate  how  our  theorem 
may  be  used  in  these  situations.  In  particular,  example  3  provides  another 
route  to  the  calculation  of  (*)  besides  that  given  in  Theorem  1  of  Glaz  and 
Naus  (1983). 


2.  Proof  of  the  Theorem. 

The  proof  uses  Che  notation  of  divided  differences  which  is  presented 
in  many  texts  on  numerical  analysis  or  the  calculus  of  finite  differences. 
Two  good  sources  are  the  texts  by  Gel* fond  and  Milne-Thomson. 

For  any  function  f  define 


f [x1,x2] 


f(x2)-f(x1) 

x2~xl 


For  n  >  2  the  divided  difference  is  defined  inductively  by 


f[x1,x2,...,xn] 


f[x2»X3 . \l-flxVx2'-“'Xn-l] 

VX1 


We  shall  need  the  following  facts.  When  x. ,x_, . . . ,x  are  all 

1  /  n 

distinct  we  may  write 


(1) 


f[x  x  ...,x  ]  -  l 
i=l 


f  f(xi>  1 


n  (x.-x.) 


3H 


i  n 


:(n-l) 


If  f  has  a  continuous  (n-1) th  derivative  f  ,  then  for  any 

x..,x0 . x  we  have 

12  n 


(2) 


(n-l)!f[x1,x2,...,xn]  ■  Ef^n_1^(  l  xiSi) 


Formula  (2)  is  an  immediate  consequence  of  a  result  due  to  Hermite  which 
may  be  found  in  section  1.6  of  Milne-Thomson. 


We  shall  abbreviate  the  notation  for  divided  differences.  For 


any  vector  x  *  (x-,x„,...,x  )'  we  let  f[x]  *  f [x, ,x,, . . . ,x  ].  Taking 
l  z  n  l  c  n 

f(u)  *  eu  (written  as  exp)  in  formula  (2)  gives  an  expression  for  the 
Laplace  transform  of  S  which  is 

(3)  E  exp(8’S)  «  (n-l)!exp[9] 

where  9  =  (8,,0, . 0  )'.  The  Laplace  transform  of  TS  may  be  obtained 

i  i  n 

from  (3)  by  inspection  and  is 


(4) 


E  exp(0 ' (TS))  =  (n-l)!exp[r'9] 


where  9  =  (01,02,. . . ,9fc) . 

Our  theorem  follows  directly  from  (4)  and  the  next  lemma  on  divided 
differences.  For  any  vector  x  *  (x^,x2, . . . ,xn> '  and  any  real  number 
u,  define  u  to  be  the  vector  obtained  by  replacing  the  i— 
component  of  x  by  u. 

Basic  Lemma:  Let  f  by  any  function  having  n-1  continuous  deriva¬ 
tives.  Suppose  d  *  (d^,d2> . . . sdn) '  satisfies  e'd  *  1.  For  all 
x  eEn  we  have 


fix] 


i-1 


V'-i,.1 


where  u  ■  x’d. 

When  uy x^ 9 x^ » • . .x^  are  all  distincty  the  lemma  may  be  verified  by 
routine  algebra  after  expanding  all  the  divided  differences  using  formula 


¥  -T  W'J  w  m  ^  m  ;  -m  '■  I'  ^  «■  -  7; 


(1).  From  (2)  it  is  clear  that  f  [x^x^ . . .  ,xn]  is  a  continuous  function 
of  its  arguments  so  that  the  lemma  must  hold  for  all  xi»x2» *  * • >xn* 

t 

Let  C  “  Td  and  u  =  0'£.  Noting  that  ( r * 0 )  =  (T.  j.)  8  we 

i >  u  s 

conclude  from  the  lemma  that 

n 

exp [ T ' 0 ]  =  l  d  exp[(r  )'0]  . 
i=l  1  1,4 

Thus  we  can  express  the  Laplace  transform  in  (4)  as  a  linear  combination 
of  similar  Laplace  transforms.  This  proves  the  theorem. 


3.  Examples . 

Example  1:  When  T  is  1  *  n  we  can  explicitly  calculate  the  distribu¬ 
tion  of  TS.  Choose  k£n  and  let  a1,a2,...,ak  be  any  real  numbers. 

We  shall  apply  the  theorem  with  T  *  (ai»a2 . an)  where  a±  “  0  for 

i  >  k. 

Assume  a^  ^  a2<  Define 


d 


1 


-  and  d„ 

a2'al  2 


al  ^ 

ara2 


Clearly  d.j+d2  *  1  and  ^lal+<^2a2  =  T^us  =  ^f^^l  0^+^2IP^r2,0^ ’ 

For  any  t  >  0,  this  implies  that 

P{[  a.S.  >  t}  *  d.Pf  l  a.S.  >  t)  +  d_P{  l  a.S  >  t> 
i  1 1  1  i*i  1  1  1^2  1  1 

where  the  sums  are  taken  over  1  <_  i  k  except  as  indicated.  Fix 
t  >  0  and  define 

6 


Q(al»a2» '  ■  *  ’  ak^  >  t }  • 


The  previous  formula  may  now  be  written  as 


(5)  Q(a^,a2, • *  * ,ak^ 


a 2*5 (a2,a3*  *  *  *  * —  al^ ^a^,a3*  *  •  • » a^) 


a2“al 


The  number  of  arguments  in  Q  is  permitted  to  vary.  Let  (z)  +  denote 
the  positive  part,  (z)+  3  niax(z,0).  It  is  easy  to  show  that 


(6) 


Q(a)  =  P{aS.  >  t}  =  (1--) 

1  E 


n-1 


Thus  when  a^ja^j.-.ja^  are  all  distinct,  repeated  application  of  (5) 
leads  to  an  explicit  formula  for  Q(a^,a2> . . . ,a^) .  This  formula  can  be 
written  as  a  divided  difference.  To  see  this  we  need  the  next  fact. 

(7)  Lemma:  Suppose  s,t,x^,x2, . . . ,xm  are  distinct  real  numbers. 
We  use  the  symbol  x  to  stand  for  the  list  x^.x^ , . . . ,xm.  Let  the 
functions  f  and  g  satisfy  g(x)  =  xf(x)  for  all  x.  Then 

sf [s,xl-tf [t,x] 
g[s,t,x]  -  - - —  . 


The  proof  of  (7)  follows  by  elementary  algebra  after  expanding  the 
divided  differences  using  formula  (1). 

From  (5)  and  (6)  we  obtain 


n-1  n-1 

a2(1‘f)  -ai(I-7-) 

z  a2  +  1  al  + 

Q(ara2)  =  - 1-^ - "  ffal’a2] 


7 


where  f(a)  =  a(l  — -)  .  Using  (5)  and  (7),  it  now  follows  by  induction 


Q(a1,a2> . .. ,a^)  g[a^,a2> . . • >a^] 


k-1  t 

with  g(a)  *  a  (1 - )  .  This  formula  has  been  demonstrated  for 

cl  T 

k  <_  n  when  the  values  a^,a2>...,a^  are  distinct,  but  it  is  easily 

shown  to  hold  for  all  values  a^,a2 . except  in  the  single  case 

when  k=n  and  a^  =  a2  =  •  •  •  =  an  =  t.  The  divided  difference  is  not 
defined  for  this  case. 

Numerical  analysis  texts  contain  prescriptions  for  evaluating 
g[a^,a2> . . . ,a^]  when  some  of  the  arguments  coincide.  For  instance, 
when  a^  =  a2=*--  =  a^  =  a  and  g  has  k-1  continuous  derivatives 
(denoted  g^)f  then 


.(k-1), 


g[a1,a2>' •* »ak]  *  (k-1)! 


This  is  a  standard  result  which  follows  immediately  from  formula  (2). 

To  use  these  prescriptions  we  must  be  able  to  evaluate  the  deriva¬ 
tives  of  g.  It  is  not  difficult  to  verify  by  induction  that 


-r(;r-)rtxm(b--)P] 

r !  dx  x 


1=0 


and  in  particular  obtain 


ijm  -d  ra  .1  _  p-i 

4(A)  [xm(b-A  i  =  y  (P)(A  (b-A 


Using  (8),  (9)  and  (11)  we  immediately  obtain 


p{  i  s  >ti  -  i  (n:1)ti(i-t)n-1~i . 

i=l  1  i=0  1 

This  formula  can  also  be  obtained  by  a  simple  direct  argument  which, 

however,  does  not  apply  in  more  complicated  situations. 

Note:  Let  be  i.i.d.  exponential  random  variables 

—  z 

with  density  e  .  If  we  let  n  -*■  00  holding  k  fixed,  it  is  easy  to 
see  that  (nS^,nS2, . . . .nS^)  converges  in  distribution  to  (Z^.Z^, . • . ,Z^) 
Therefore,  substituting  t  =  x/n  in  formula  (8)  and  letting  n  -*•  » 
yields 

k 

(12)  H  l  ajLZi>x}  =  f[a1,a2 . ak] 

,  , ,  .  k-1  -x/a 

where  f(a)  »  a  e 

Partial  Divided  Differences. 

In  order  to  generalize  the  result  of  example  1  we  need  to  introduce 

divided  differences  for  functions  with  more  than  one  argument.  The 

definitions  will  be  presented  in  terms  of  a  function  f  having  two 

arguments,  f  =*  f(x|y).  The  definitions  are  similar  for  functions  with 

more  than  two  arguments.  Let  u  and  w  represent  arbitrary  lists  of 

real  numbers  uisu  , ...  u  and  w, ,w„ , . . . ,w  .  Define  the  divided 
12  p  12  q 

differences  recursively  by 

f(x2|y)-f (x1|y) 


f[x1,x2|y] 


9 


f[x|y1,y2l  = 


f(x|y2)-f(x|yi) 


y2“yl 


f[x1,x2,u|w]  * 


f[x2,ujw]-f [x1,u|w] 
X2“X1 


and 

f [u|y2»w]-f [u|yi,w] 


It  may  be  shown  that  the  divided  differences  do  not  depend  on  the  order 
in  which  the  differences  are  taken.  In  particular,  for  each  fixed 
value  of  x  we  may  define  the  function  f^  by  fx(y)  =  f(x|y)  and 
write 


f [u|y1,y2,w] 


(13)  f [x^,x2, . , . ,Xp| y^,y2> . . . ,y  ] 

=  g[x1,x2, . . . ,X  ]  where 

g(x>  =  fx[y1>y2’'--’yq]  • 

w.  the  values  x^,x2>...,Xp  are  distinct  and  y^.y^, . . . .y^  ari 
also  d  inct,  the  divided  difference  is  given  by  the  formula 

(14)  f[xlt...,xp|ylt...,y  ] 

f(xi|yi> 

(x  -x  )  n  (y  -y  ) 

1  k  3  1 


i=l  j=l 


n 

k*i 
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1  1  0  0  0  0 
0  110  0  0 
0  0  1111 


1  0  0  0  0  0]  M  i1  0  0  0  0) 

0  1  1  1  1  lj  \o  lj  +  (o  1  1  1  1J  ‘ 


Putting  this  all  together  gives 


[l  0  0  0  0  o|  [l  0  0  0  o]  [l  0  oj  [i  o' 
{A}  =  2  \o  1  1  1  1  lj  +  2  (o  1  1  1  lj  '  2(o  1  1/  '  (o  lj  • 


Evaluating  these  terms  using  formula  (22)  leads  to  the  final  answer 


(l-2t )^_1+2 (n-l)t (l-2t)+~2+4 (n"!) t2 (l-2t)^'3+4 (n"!) t3 (l-2t)^-4 


+  2(n~1)t4(l-2t)J"5. 


Tables  given  by  Naus  (1966)  provide  a  numerical  check  for  this 

answer.  Naus'  notation  differs  from  ours.  Using  our  formula  (with 

n*8)  to  calculate  l-{A}t  for  the  values  t  =  .10,  .20,  .30,...  leads 

to  the  numbers  given  in  Naus'  Table  2  for  the  case  N=7,  n=4. 

Note:  As  in  the  note  after  example  1,  we  can  transform  our  results 

about  spacings  into  results  about  i.i.d.  exponentials  by  substituting 

t  *  x/n  and  letting  n  -*•  °°.  Let  Z^>Z2»*‘*»Z6  be  with  density 

-2 

e  .  The  result  of  example  6  implies  that 


P{Z^+Z2+Z^  >  x,Z2+Zyt-Z^  >  x,Z|j+Z^+Z^  >  x,Z^+Z^+Z^  >  x} 


e~2x(l+2x+4(x2/2! ) +  4(x3/3!) +  2(x4/4! ) } 


22 


evaluated  using  formulas  (20)  or  (21).  However,  it  is  convenient  to 
further  simplify  some  of  the  terms. 


1  0  0  0  o' 

1  o  o  o' 

i  o  o  o' 

,01000 

►  =  2  . 

0  10  0 

►  -  • 

0  10  0. 

0  0  2  1  1 

0  0  2  1 

0  0  1  1 

and 


1  0  0  o' 

1  0  o' 

1  0  o' 

.0100 

0  0  2  lj 

’  2  1 

0  10 

0  0  2 

0  10. 

0  0  1, 

because  in  both  cases  2y^-y^  *  (0,0,0)'.  Putting  all  of  these  results 
back  in  our  original  expression  gives 


{A} 


1  0  0 

1  0  o' 

1  0  0  o' 

16  , 

0  10 

f  -  10  . 

0  10 

r  - 5  ’ 

0  1  0  ol 

[0  0  2, 

0  0  1 

.0011 

Using  formulas  (20)  or  (22)  to  evaluate  these  terms,  we  obtain 


{A}  -  16(l-(5/2)t)“_1-  15(l-3t)"_1-  5(n-l)t(l-3t)"-2  . 


Example  6: 

Our  final  example  will  be  to  compute  {A}t  for 


A 


1  1  1  0  0  o' 

0  i  1  1  0  0 
<  ► 
0  0  1110 

0  0  0  1  1  1 
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Y^  will  be  used  generically  to  denote  the  i — -  column  of  a  matrix. 
We  begin  by  noting  that 

l-l+l-l+l  -  1  and 
Y  l-^4’*  3~Y4+Y  5  =  <0.0.0, 0.2)' 

so  that  our  theorem  expresses  {A}  as  an  alternating  sum  of  five  terms, 
successively  replacing  each  column  by  (0,0, 0,0, 2) ' .  By  using  properties 
(a)  and  (d)  to  simplify  these  terms  and  then  combining  terms  which  are 
seen  to  be  equal  upon  using  properties  (b)  and  (c),  we  obtain 


Noting  that  the  columns  of  the  first  term  satisfy  *  (0,0, 0,1)', 

the  first  term  can  be  rewritten  as 

0  10  0  0 
0  0  110 
0  0  0  1  1 
2  0  0  0  1 

where  again  we  have  simplified  and  rearranged  the  matrices  using  (a) 
through  (d).  All  the  terms  in  our  expression  for  ( A }  can  now  be 


'l  0  0  0  Oj  fl  0  Oj  fl  0  0  0 

Jo  100  oj*  -  <  o  1  ol  +  <  o  1  00 

0  0  2  1  lj  [o  0  l)  [o  0  2  1 
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Using  (b)  we  can  rewrite  this  as 


[0  1  0)  fl  0]  (l  0  0 

-  -  •  + 

U  o  lj  i|  to  i  i] 


Jl  0  0]  _  fl  0)  . 

fo  1  lj  lo  lj 


These  terms  can  now  be  evaluated  using  formulas  (20)  and  (22)  to  yield 


2[(l-2t)"-1  +  (n-l)t(l-2t)"'2]  -  ( 1— 2 t ) 


n-1 


(l-2t)"-1  +  2(n-l)t(l-2t)"-2 


Example  5: 

A  modification  of  the  previous  example  is  to  compute  {A}t  with 


110  0  0 
0  110  0 
^00110^ 
0  0  0  1  1 
1  0  0  0  1 


Changing  the  dimension  of  A  from  A  to  5  produces  a  considerable  change 
in  the  complexity  of  the  calculation  and  the  form  of  the  answer.  When 
n-5,  the  probability  (A}t  has  an  interpretation  similar  to  that  in 
the  previous  example. 


The  distribution  of  the  scan  statistic  (and  related  matters)  has  been 
extensively  investigated  by  Naus  and  his  co-workers  Wallenstein, 
Huntington  and  Glaz. 


We  shall  now  compute 


P{Sj+S2  >  t,S2+S2  >  t.S^+S^  >  t.S^+S^  >  t}  =  (A)t  with 


110  0 
0  110 
0  0  11 
10  0  1 


When  n»4,  can  be  viewed  as  the  spacings  between  four 

points  placed  uniformly  on  a  circle  with  circumference  equal  to  one 
and  thus  {A}fc  is  the  probability  that  no  arc  of  length  t  contains 
more  than  two  of  the  four  points. 

Let  denote  the  i—  column  of  A.  Since  Y^_Y2+Y3  “  ^4  anc* 

1-1+1  *  1,  our  theorem  implies  that 


{A}  *  {A.  }  -  {A„  ,  }  +  (A.  w  } 


0  10  0 
0  110 
10  11 
10  0  1 


10  0  0 


0  111 
110  1 


1  1  0  01 


>0  0  1  ol  +  io  1  0  oj 


0  0  11 
10  11 


By  applying  property  (a)  and  then  property  (d)  to  each  term  we  obtain 
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Some  Notation  and  Elementary  Properties. 

Let  A  be  any  matrix  having  at  most  n  columns.  Take  r  to  be 
the  matrix  with  n  columns  obtained  by  padding  A  with  columns  of 
zeros;  r  *  (A|0).  Let  k  be  the  number  of  rows  in  A  and  define 
X  ■  (X^,X2» • • • ,X^) '  by  X  *  TS.  For  any  value  of  t  we  define 

{A}'t  =  PfX^  >  t  for  all  i}  . 

When  the  value  of  t  is  held  fixed  in  an  argument,  we  shall  delete 

the  subscript  and  just  write  {A}. 

Our  theorem  and  its  consequences  can  sometimes  be  used  to  evaluate 

{A}.  The  following  elementary  properties  are  useful  in  the  manipulations, 
til 

(a)  If  the  i—  row  of  A  dominates  (is  componentwise  greater  than 

t  ll  t  ti 

or  equal  to)  the  j —  row,  the  i—  row  can  be  deleted  without 
changing  the  value  of  {A}. 

The  value  of  {A}  does  not  change  when 

(b)  the  columns  are  permuted, 

(c)  the  rows  are  permuted,  or 

(d)  a  column  of  zeros  is  deleted. 

The  next  three  examples  are  calculations  which  in  some  simple  cases 
(and  with  the  appropriate  choice  of  n)  give  the  probabilities  of  certain 
types  of  clustering  on  a  circle  or  interval.  These  sorts  of  calculations 
may  be  used  to  study  the  distribution  of  the  scan  statistic,  a  test  for 
the  presence  of  non-random  clustering  (see  Naus  (1966b)  or  Cressie  (1977)). 
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The  quantity  (  )  is  the  usual  multinomial  coefficient. 

i *  r2 ' '  ”  ’  p 

The  expression  (22)  is  easily  derived  by  direct  methods  which,  however, 
do  not  apply  in  more  complicated  situations. 

Example  3: 

We  now  show  how  to  compute  P(D  when  T  =*  (y^)  is  2*n  and 
*  0  or  1  for  all  i  and  j .  Let  y^  denote  the  i—  column  of 
T.  Define  r,  s  and  t  to  be  the  number  of  columns  in  T  which  are 
equal  to  (1,0)',  (0,1)'  and  (1,1)'  respectively. 

Since  P(T)  is  invariant  under  permutation  of  the  columns  of  r, 
P(F)  depends  only  on  r,  s  and  t  so  that  we  may  define 
P(r,s,t)  -  P(D .  Assume  r  >  0,  s  >  0  and  t  >  0,  then  without  loss 
of  generality  we  may  take  yx  -  (1,0)',  y2  -  (0,1)'  and  y,  *  (1,1)'. 
1+1-1  *  1  and  YJ+Y2-Y3  *  (0,0)'  *  0,  our  theorem  gives 

*<n  -*<rli?)+  P(r2>0)-P(r,  0)  . 

This  may  be  rewritten  as 

(23)  P(r,s,t)  -  P(r-l,s,t) +P(r,s-l,t)  -  P(r,s,t-1)  . 

Repeated  application  of  this  recursion  allows  us  to  express  P(r,s,t) 
as  a  sum  of  distributions  P(i,j,k)  in  which  i-0  or  j  *  0  or 
k>*0.  These  distributions  are  easily  handled  by  other  methods. 

A  generalization  of  (23)  was  obtained  by  a  different  method  in  Huffer 
(1982). 


Since 
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p{Zi>ti  for  all  i}  =  C[x1(x2|  ••  •  |x  ] 


P  k  -1 


(21)  where  £  (yL  jy2  |  •  •  •  |  y  )  =  (Eh<y1S1 1  y2s2 1  *  -  - 1  y  S  ) )  n  y^ 


- 1  wr1  s  n1'1  • 

i-l  +  i-l  1 


When  the  components  of  the  vectors  are  distinct,  the  divided 

differences  can  be  evaluated  by  an  obvious  generalization  of  (14).  When 
the  components  of  each  vector  are  all  equal,  we  can  calculate  the  divided 
differences  by  generalizing  formula  (15).  For  1  ^  i  <  p,  let 
e^  «  (1,1,..., 1)  have  ki  components.  It  is  easy  to  show  that 


ctyi;ily2e2|.-.|yp!pl  -  opiTT J ^ Cy 1 1 y2l - •  •  I yp) 


With  £  as  given  in  (21) ,  these  partial  derivatives  are  obtained  by 


repeated  use  of  (11) .  Setting  y^  =  y2 ; 


yp  »  1  in  the  resulting 


expression  gives  a  formula  for  the  joint  distribution  of  nonoverlapping 
sums  of  uniform  spacings. 


JL 

P{  £  S  >  t^  for  1  <  i  <  p} 


j^l-H 


I  (  n_1 

r  L  \  r, , r- , . 

rl’r2’ ♦ *  * ,r_)  I  2’ 


p  n-l-I  r.  p  r. 

r  )(i-  It.)  1  1  n  t  J 

’’  p  i-l  +  j»l  J 


°£ri£ki-i 


where  m^  -  0  and  m^  »  k^  for  1  £  j  £  p.  The  summation  is 

taken  over  all  p-tuples  of  integers  satisfying  0  <  r.  <  k.-l  for  all  i. 


This  result  can  easily  be  extended  to  handle  p  >  2.  We  shall 
state  only  the  following  special  case.  Assume  T  has  a  block  diagonal 
fora 


h  o 

\ 

o 

where  each  x.^  is  a  row  vector  having  components.  Let 


Z 


(z1,z2,...,z  r-  rs  , 


choose  real  numbers  t. ,t„,...,t  and  define 

i  i  p 


1  if  z^>  ^  for  all  i  , 
0  otherwise  , 


so  that  Eh(Z)  »  P{Z^>  for  all  i}  . 

It  is  well  known  (see  Feller  volume  II,  problem  23  of  chapter  I) 

that 


(20) 


>  t1,s2  >  t2, 


•vv 


p 

(1- 1 

i*l 


ti>r1 


Repeated  use  of  (16)  followed  by  an  application  of  (20)  leads  to  the 


formula 


» 


X1  X2 

*k<xl,x2***  *,xk*  =  <x1-x2)^k-l(xl,X3’  ‘  ’ '  ,xk)  +  (x2-x1)\-l(x2,x3’ - '  *  »*k* 


By  induction  as  in  example  1  we  now  obtain 


♦  k(xi**2»”- V  =  *fx1»x2»  •  •  •  »xk] 


with  4>(x)  =  xk  %^(x) 


Assume  now  that  p  =  2  and  F  has  a  more  specialized  form 


x  0  0 


lolyloj 


where  y  is  an  m-dimensional  row  vector,  h  is  now  a  function  of  two 
arguments,  h  =  h(xjy).  Applying  (16)  to  the  quantity  ^(x)  viewed 
as  a  function  of  y  yields 


i/^Cx)  =  0[y1,y2,..  .,ym] 


with  0(y)  ■  ym  Eh(xS^)yS2) 


Combining  (16)  and  (18)  leads  to 


Eh(rS)  =  £[x|y] 


with  C (x| y)  =  xk  1ym  1Eh(xS1[yS2) 


whenever  F  has  the  form  in  (17). 
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This  is  easily  proved  using  formula  (1)  and  the  representation  (13). 

When  some  of  the  values  coincide,  limiting  arguments  yield  prescriptions 
(in  terms  of  derivatives)  for  evaluating  the  divided  differences.  For 
example,  when  x^  =  x  for  all  i  and  y  = y  for  all  j,  then 


(15) 


f  ,  .  .  .  ,  X^  I  y  ^  j  *  *  *  » 1 


(p-l)!(q-l)!  (3x*  V  • 


This  result  follows  from  (9)  and  (13). 


Example  2: 

Let  T  by  any  p  x  n  matrix  which  can  be  written  in  the  block  form 


I 

b 

0 

0 

A 

0 

where  x  =  (x^,x2 , . . . ,x^)  is  a  k-dimensional  row  vector  and  0  denotes 
a  matrix  or  vector  of  zeros.  For  an  arbitrary  bounded  function  h  on 

p 

R  we  define 


'l'k(Xi,x2,...,xk)  =  Eh(TS)  . 

At  first  we  shall  consider  b  and  A  to  be  fixed.  The  subscript  k 
is  included  in  the  notation  because  we  shall  be  varying  the  dimension 
of  the  vector  x. 

Assume  that  x^,x2>...,xk  are  all  distinct.  Using  our  theorem 


exactly  as  in  example  1  we  obtain 
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